clear all
load ParEst_eh1_time0_Mex_guys ParEst_eh1_time0_Mex_guys
%load ParEst_eh1_time0_Mex ParEst_eh1_time0_Mex
load trans_time1_age0_eh1_comp1.out
Chute =ParEst_eh1_time0_Mex_guys;
% Chute=[ParEst_eh1_time0_Mex(1), ParEst_eh1_time0_Mex(2), ParEst_eh1_time0_Mex(3), ParEst_eh1_time0_Mex(4), ParEst_eh1_time0_Mex(5), ParEst_eh1_time0_Mex(6),...
% ParEst_eh1_time0_Mex(17),ParEst_eh1_time0_Mex(18),ParEst_eh1_time0_Mex(19),ParEst_eh1_time0_Mex(20),ParEst_eh1_time0_Mex(21),ParEst_eh1_time0_Mex(22),...
% ParEst_eh1_time0_Mex(23)];
HandleLo = @Principal_A2_n_1;

D_mom=trans_time1_age0_eh1_comp1; % Read transition moments
Dij=D_mom';

Dnf_1	 =	Dij(1);
Dni_1	 =	Dij(2);
Dfn_1	 =	Dij(3);
Dff_1	 =	0;
Dfi_1	 =	Dij(4);
Din_1	 =	Dij(5);
Dii_1	 =	0;
Dif_1  =	Dij(6);

Lower =[Dfn_1,Din_1,Dnf_1,Dni_1,Dfi_1,Dif_1];
LoBo =zeros(1,13);
LoBo = [Lower, 0.1, 0.1, 0.1, 0.1, -10, -10, -10];
HiBo = [.81 .81 .81 .81 .81 .81 10 10 10 10 -4 5 0];
                 options = saoptimset( ...
                     'PlotFcn',{@saplotbestf,@saplottemperature,@saplotf,@saplotbestx,@saplotx},...
                     'TolFun',5000,'ReannealInterval',6, ...
                     'InitialTemperature', 14*10^2*[1 1 1 1 1 1 10 10 10 10 6 15 10], ...
                     'TemperatureFcn', @temperaturefast, 'StallIterLimit', 1000); %default 14*10^2
[x,fval,exitflag,output]  = simulannealbnd(HandleLo, Chute, LoBo, HiBo, options)
ParEst_eh1_time0_Mex_guys = x;
save ParEst_eh1_time0_Mex_guys ParEst_eh1_time0_Mex_guys
